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Abstract 

The resistive switching phenomenon in MgO-based tunnel junctions is attributed to the effect 
of charged defects inside the barrier. The presence of electron traps in the MgO barrier, that can 
be filled and emptied, locally modifies the conductance of the barrier and leads to the resistive 
switching effects. A double-well model for trapped electrons in MgO is introduced to theoretically 
describe this phenomenon. Including the statistical distribution of potential barrier heights for 
these traps leads to a power-law dependence of the resistance as a function of time, under a 
constant bias voltage. This model also predicts a power-law relation of the hysteresis as a function 
of the voltage sweep frequency. Experimental transport results strongly support this model and 
in particular confirm the expected power laws dependencies of resistance. They moreover indicate 
that the exponent of these power laws varies with temperature as theoretically predicted. 
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I. INTRODUCTION 



Resistive switching effects^ have been studied since the 70s in range of insulating oxides 
such as Ti02 or AI2O3 for instance^^"-. This interest has been renewed for a few years, 
since giant and reproducible effects were observed in perovskites^ as for example SrTiOs 
doped with chromium. It makes these materials good candidates for a new generation of 
memories. Yet, the underlying physical mechanisms are still unclear and different hypotheses 
have been put forward. Electro-migration of dopants or oxygen vacancies along filaments 
could reversibly create conducting path across the insulating layer-i^. Another hypothesis^ 
suggests the accumulation of charges at the electrode/insulator interface, which depends on 
the applied bias, and thus changing the Schottky barrier height. 

The switching effect is in most cases studied on relatively thick films, on the order of 
100 nm thick, but it has also been observed in some systems with a thin barrier allowing 
tunnel transport^i"— . We showed for instance that MgO tunnel barriers with a few atomic 
layers of chromiunt^^ or vanadiunt^^ at the MgO interface exhibited reproducible switching 
effects in Fe/Cr/MgO/Fe or Fe/V/MgO/Fe systems. This was attributed to the creation of 
oxygen vacancies in MgO at the interface with these "dusting" layers. These defects locally 
open extra channels in parallel with the "standard" tunnel transport through the MgO 
barrier. Moreover, these systems exhibit an interesting behavior, with a relaxation of the 
conductance on long time scales -on the order of minutes. We indeed observed logarithmic 
relaxation of the conductance under a constant bias voltage U. A strong influence of the 
voltage sweep frequency on the hysteresis in I{U) curves was also observed. 

In a recent article^^ we proposed a phenomenological model accounting for the relaxation 
of the conductance with time and the effect of sweep frequency on the hysteresis in I{U) 
curves. This model did not discriminate between both hypotheses: either electro-migration 
of oxygen vacancies could create local conducting paths in the barrier, or the accumulation of 
charges inside the barrier could modify its potential height and thus the tunneling transport. 
Within the scope of this general model, both hypotheses could lead to the same mathematical 
expressions. In this model, the relaxation of conductance with time was expected to be 
exponential and not logarithmic. 

The present article partly justifies this phenomenological model and supports the hy- 
pothesis of electron trapping as the microscopic origin of the modification to the tunneling 
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conductance. We indeed observe telegraphic noise in the conductance which is interpreted 
as the sign of electron trapping and untrapping due to its low activation energy. We then 
propose a double-well model to account for these trapping events. In this new model, double 
wells are characterized by a random barrier height between the two trap states. Assuming 
an exponential tail for the distribution of these barrier heights, we obtain a power law de- 
pendence of the current hysteresis AJ as a function of the voltage sweep frequency of I{U) 
curves. Moreover, the exponent of the power law obtained in this model is proportional 
to the temperature. We performed systematic measurements as a function of temperature 
which clearly confirm these theoretical predictions. 

II. EVIDENCE OF ELECTRON TRAPPING 
A. Telegraphic noise 
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FIG. 1: (color online) I{U) curve measured at 80K showing hysteresis above +170 mV and a low 
current under a negative applied voltage. 

We present results obtained on a Fe(20nm)/V(1.2nm)/MgO(3nm)/Fe(5nm)/Co(15nm) 
sample grown by molecular beam epitaxy. Details of the growth are given in Ref. [15] to- 
gether with the details of the micron-sized junction processing. Electrical measurements are 
performed with a conventional four-point DC technique. The reference of positive voltage 
was taken as the top electrode (with no vanadium). I{U) curves on such samples were 
already shown at room temperature in Ref. [15] and exhibited systematic hysteresis. At low 
temperature, the I{U) curves are still hysteretic as can be seen in Fig. 1. We have to note 
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a threshold (in the order of + 170 mV) below which the hysteresis is absent. Moreover, 
the junctions behave as a rectifier at low temperature: the current under negative bias be- 
comes much smaller than under positive bias. This might be attributed to the asymmetric 
potential barrier in presence of the vanadium layer. Indeed, the tunnel transport through 
monocrystalline MgO barrier is dominated by electrons having the Al symmetry^^ for which 
vanadium represents a large potential barrier— - more than 4.2 eV- due to its band struc- 
ture. It thus leads to an asymmetric barrier with rectifying characteristics. This point has 
to be further studied. In the following, we will just show results obtained at a positive bias. 
We observe for a low constant bias voltage (less than the threshold value of +170 mV), a 
telegraphic noise (see Fig. 2(a)), proving a bistable conductance of the junction. Provided 
that the voltage is lower than the threshold, almost no relaxation of the conductance value 
is observed on long time scales: the average level of conductance remains almost constant. 
We have to stress the fact that Fe/MgO/Fe samples which do not show resistive switching 
effects or relaxation with time^^ do not show this telegraphic noise either. This supports the 
idea of a correlation between this observed noise and resistive switching mechanisms. 



B. Energy levels of electronic traps 

In the case of samples showing telegraphic noise, by slightly changing the applied dc 
bias U, we modify the occupancy rates between the low and high resistance states labeled 
respectively (1) and (2). This enables us to plot the ratio of occupancy states ^ as a function 
of U and to fit this ratio as an exponential dependence on U (see the linearity on the Log 
plot of Fig. 2(b)). Supposing that the time occupancy for both states follows an Arrhenius 
law^^i^ we obtain that 

-1 = Ke^, (1) 

Tl 

where AE = i^i — £^2, K is a constant, ks is the Boltzmann constant and T is the tem- 
perature. The energy of both states can moreover be written (see Fig. 4): Ei = + aU 
and E2 = E2 — aU where U is the applied voltage and Ef is the energy of state i in the 
absence of applied voltage. Thus, from the slope of kBT\\i{—) we can extract ^ = 2a 
which corresponds to the voltage-dependent part of the energy difference between the two 
states. This yields = 135 meV/Volt for the measurement made at 80K. We found 
similar values on other junctions of the same sample. 
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By extrapolating our plot of Fig. 2 to [/ = 0, we can have access to = E2 — 
provided we make the hypothesis that K = 1. We then find that E^ = —19 meV. 

This value and its voltage dependant part a are very low compared to reported values 
obtained thanks to telegraphic noise in other devices^^: for instance Eq = 250'meV in Cu 
doped SeGe resisting switching systems^^. It is all the more striking as our values are 
obtained in MgO in which atomic displacements require high activation energies^Si^in the 
order of 2 eV. This is two orders of magnitude higher than what we observe. It thus supports 
the hypothesis of charge trapping instead of atomic displacement in order to explain the 
resistive switching in our system. 

It is tempting to make the same type of measurements as a function of temperature. 
Unfortunately, obtaining such a dataset as a function of T was not possible: in some cases, 
no telegraphic noise is observed after increasing the temperature and stabilizing it -which 
takes more than 15 minutes- . Moreover, we can still observe such telegraphic noise for the 
new temperature, but associated to other metastable states. This problem is not so crucial 
when changing the voltage at a given temperature: this is made rapidly, by small voltage 
steps, thus "following" the two states. Notice that the behavior shown on figure 2(b) is an 
indication that the same two states are considered when sweeping the voltage. This behavior 
was not observed as a function of temperature, suggesting that the telegraphic noise was 
associated to different states when modifying the temperature. 

This voltage dependence suggests that charges can be trapped in the MgO barrier or at 
the interface: whether they are trapped or untrapped, the potential height of the tunnel 
barrier might be modified, thus infiuencing the probability for electrons to tunnel from one 
electrode to another. A similar phenomenon has already been observed in MOSFETs^^ 
channels below the grid insulating oxide, or in thin Josephson^S junctions. In both cases, 
oxygen vacancies in the oxide create charged defects whose charge fiuctuates over time, 
leading to telegraphic noise in electronic transport as observed in our case. 

In our case we can suppose that the electric charging of the barrier locally modifies its 
potential height due to electrostatic effects and thus changes the tunnelling probability of 
electrons close to this trap, as in the case of Coulomb Blockade. This effect might lead to 
strong conductance changes when the trap is located on a hotspot: it was indeed shown^Si"— 
that the tunnel transport through such thin barriers is not homogenous but dominated by 
some hotspots. They can for instance be due to a locally thinner barrier because of the 
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FIG. 2: (color online)(a) Resistance measurement at 80K under a constant voltage U=+120mV 
showing telegraphic noise. This is characteristic of a bi-stable single defect oscillating between 
two states with two different conductance values. Inset: Histogram of the curve showing both 
populations on states 1 and 2 with two gaussian fits, (b) (squares) Experimental ratio of residence 
times as a function of voltage, calculated from histograms. T2 corresponds to the high resistance 
state, (full line) Linear fit. 

roughness of the oxide. This rare events effect explains that a single electronic trapping 
could lead to 2 % changes of the junction conductance, provided it is located on such a spot. 

Moreover, the value obtained for can be linked with the position of the defect inside 
the MgO barrier. Let us indeed suppose that the involved charged particle is an electron and 
that one trap position is at the V/MgO interface and the other at the MgO/Fe interface. 
The change of potential for the electron moving from one trap to another would he q ■ U 
where q is the electron charge and U the applied voltage. The value of should then 
be leV/Volt which is higher than what we observe. This proves that the two traps are not 
located at the junction interfaces of the MgO barrier, but rather at a distance d of each 
other with 

dAE 

~ dU 
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where dugo is the MgO barrier thickness. We find d = 0.4nm. This would be consistent with 
the creation of oxygen vacancies at the lower interface close to the vanadium layer: 0.4nm 
gives the order of magnitude of the thickness of the faulted MgO layer containing traps. 

Looking at Fig. 2 we observe that a higher positive voltage favors the higher resistance 
state. Within our polarity convention, a positive voltage corresponds to electrons moving 
from the bottom -with vanadium- interface to the top interface. In the case of trapping it 
thus suggests that a positive voltage favors charging the traps inside MgO near the V/MgO 
interface, and leads to a decrease of the conductance across the barrier. 

Another point has to be stressed: at positive voltage the voltage-dependent part of the 
activation energy of traps and the constant part have opposite signs. Two regimes can thus 
be distinguished relative to 



which is here of 140mV. This value should be compared with the threshold value for hystere- 
sis in I{U) curves, close to 170mV. Above this value, the volt age- dependent term dominates, 
leading to a partial filling of MgO traps -and to an increase in the junction resistance. Below 
this voltage threshold value the constant term dominates, the trap states depopulate and 
the hysteresis disappears. 

III. RELAXATION WITH TIME AND ROLE OF THE VOLTAGE SWEEP FRE- 
QUENCY 

We now turn to a higher voltage regime, i.e, with U > +170mV. We showed in a previous 
article^ that the relaxation of the resistance under a constant voltage was nearly logarithmic 
at room temperature. Here we have performed measurements of this relaxation at low 
temperature -from lOK to 200K. As seen in the data of Fig. 3, the behavior remains the 
same at low temperature as shown on Fig. 3. We have to stress that no telegraphic noise is 
observed in this high voltage regime. A Log-Log plot of the conductance G as a function of 
time yields a linear plot, corresponding to a power-law with an exponent m ^ 1. This value 
strongly depends on the estimated value of Go which is the constant part of the conductance, 
corresponding to the asymptote on G(t) curves. For instance, as shown on Fig. 3 at 80K, 
m ~ 0.06 ± 0.01 if we take Go = 1.2.10~'^S' , whereas we obtain m ~ 0.16 ± 0.04 if we take 
Go = 3.6. lO"^^. Within this [0.06 — 0.16] range all fits are correct. 




(3) 
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Nevertheless, this low value of m indicates that the curves can be regarded to first ap- 
proximation as almost logarithmic, as we did in Ref. [15]. 




FIG. 3: (color online) (a) (dots) Conductance G of the junction at 80K as a function of time under 
a +400 mV bias after applying for 3 minutes a -300mV voltage, (squares) Fit corresponding to 
G = Gq + at~"^ where Go = 1.2.10~^S' is the evaluated non switching part of the conductance, 
a is a constant and m ~ 0.06. Inset: G — Gq = at~^ on a Log-Log scale, with two different 
values for Gq yielding two different exponent m as explained in the text, (b) Hysteresis in current 
A/ measured at +0.4V at 80K as a function of the bias sweep frequency. A/ scales as uj"^ with 
m' = 0.22 ± 0.05 (straight line). 

Furthermore, in the case of dynamical measurements -i.e., when making I{U) 
measurements- we have shown^ that the hysteresis observed on I{U) curves depends on 
the frequency of the voltage sweep: the A/(ci;) curve shows an increase at low frequency 
followed by a slight decrease at high frequency. If we focus on the low frequency regime 
(see Fig. 3b), A/ follows a power law as a function of at T = 80K, we find an exponent 
m' ^ 0.22 ±0.05. 
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IV. RANDOM BARRIER DOUBLE- WELL MODEL 
A. Definition of the model 




FIG. 4: (a) Electronic potential landscape for an electron on two traps, without applied voltage, 
(b) With an applied positive voltage U. 

The very slow relaxation of the conductance indicates that the system does not possess a 
single relaxation time scale, but rather a broad distribution of time scales, suggesting that 
disorder effects may play an important role. It has been known for a long time that the 
presence of disorder may strongly affect the electrical properties of materials, for instance in 
ionic conductors like "HoUandite"— , or in amorphous insulating materials like As2Se^^, 
and models based on random distributions of barrier heights have proved useful to account 
for the behavior of such systems^i^i^. In this section, we propose a simple double-well 
model with a random barrier between the two wells in order to describe our experimental 
data. 

We develop this model under the hypothesis of electron charging of traps in MgO as 
suggested by experimental observations shown above. It would nevertheless give exactly the 
same mathematical results under the hypothesis of atomic drift of atoms, locally modifying 
the conductance. The double-well can indeed correspond to two positions of the involved 
ion inside the tunnel barrier, yielding two different values of the local tunnel conductance. 

We consider that the electronic conductance G results from many independent parallel 
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conductance channels. Some of these channels are 'standard', and give altogether a contri- 
bution Go to G. The other channels are modeled by double-well potentials in which electrons 
can be trapped. These wells are assumed to be separated by an energy barrier of random 
height. In a given channel, the two potential wells have an energy El = + aU and 
E2 = -£^2° ~ CiU respectively {E^ < E2), where U denotes the electric potential and a is 
an effective electric charge. To each well is also associated a given tunneling conductance, 
denoted as gi and g2, with g2 > gi- For simplicity, we assume that E^, E2, a, U, gi and 
g2 have the same values in all the channels. In contrast, the energy barrier varies from one 
channel to another. With an appropriate choice of the energy reference, we set E^ = —Sq 
and E2 = Sq. An electron going from the first well to the second one has to cross an energy 
barrier 

A^i = W + So-aU, (4) 
which defines W (see Fig. H]). In the opposite direction, the energy barrier is 

AE2 = W-So + aU. (5) 

Hence W can be interpreted as the average barrier between the two wells. We consider W 
as a random variable, and denote as p{W) its probability distribution. 

Assuming that a number Uc of non-standard channels are present in the system, the total 
conductance at time t is given by 

G{t) = Go + ricplit) gi + n,{l - p{{t)) g2, (6) 

where pi{t) denotes the average occupancy rate of the well of conductance gi (the average 
being performed on the different channels, or equivalently, on the statistics of the barrier 
W). At a temperature T, the mean time to cross the barrier AEj is given by an Arrhenius 
law 

r, = roe^^^/'=-^, (7) 

where tq is a microscopic time characterizing the vibrations at the bottom of the wells. We 
introduce the occupancy rate pi{t,W) of the first well, given the barrier W. The average 
occupancy rate is then obtained by averaging over the barrier W, namely 

Wit) = {piit,W))w. (8) 
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The evolution equation for pi(t, W) reads 



^{t, W) = --p,it, W) + 1(1 - p^{t, W)), (9) 

which can be rewritten, using Eqs. @ and ([5]) as 

^(t,W) = le-^/^-^ (10) 
at To 

X [e^WAfl^ _ 2pi(t, W) cosh{S {t)/kBT)] , 

with S{t) = So-aU{t). 

B. Response to an electric potential step 

The relaxation of pi(t,W) after a step in the electric potential U(t) = UoQ{t) (where 
6(t) is the Heaviside function) is readily calculated, yielding 

Pi{t, W) =pf + A exp (-7t e-^/'^^^) , (11) 

with , A and 7 given by 

^1 ^ 1 + e-2(So-aC/o)/fcsT (-'-^) 

^ ^ 1 + g-25o/fcBT ~ 1 _|_ g-2{So-a(7o)/fcsT 

7 = - cosh , ^ . (14) 



To V ksT 
Averaging over the barrier W yields 

p-i{t) = pt + A(exp(-7t 6-^/'=-^))^/. (15) 
The average of the exponential term reads 

/•oo 

(exp(-7te-^/'^^))iy = / dWp{W) (16) 

xexp(-7te-^/'=^^) 

where Wmm is the minimum value of the barrier W. In order to compute explicitly this last 
average, we need to choose a specific form for p{W). Following the standard literature on 
trap and barrier models^iiSSj we consider a distribution p{W) with an exponential tail, 

p(W^) ~ Ce-^/^o, W^oo, (17) 
11 



where C > is a constant. Such a form can be justified for instance on the basis of extreme 
value statistics^. If p{W) is purely exponential, C is given by 

C = Wo-^e^-'"/^°. (18) 

Making the change of variable z = ^-"^/^bT ^ obtain for large time t 

{eM-lte-^^'^'))w - ^^^^ (19) 

with = ksT /Wq, and where r(a;) = dyy^~^e~^ is the Euler Gamma function. Accord- 
ingly, we have 

From Eq. ([6]), we thus find that the conductivity G{t) relaxes as a power law to its 
asymptotic value, with an exponent fi proportional to the temperature. If the temperature is 
small, namely /i ^ 1, then the relaxation is approximately logarithmic over a significant time 
window. We underline that, within our model, the evolution of resistance with time is driven 
by the populations of states 1 and 2 which varies continuously with time. The fiuctuation 
of resistance with time being averaged on a large number of defects, no telegraphic noise is 
expected. 



C. Response to a periodic excitation 

We now turn to the case of a small periodic excitation U{t) = uocos{ut), such that 
auo ^ kBT. We first consider a single channel, with a fixed barrier W. Starting from 
Eq. (fTOj) , we look for a solution of the form 

Pi{t,W)=pU^pl{t,W) (21) 

and we linearize Eq. ffTOj) with respect to the small parameter auQ/ksT. The zeroth order 
equation yields 

^1 = l + e-25o/fcsT- (22) 

At first order in auo/ksT, we get 

r^e^/^BTM. = ^2p\{t, W) cosh ^ (23) 



— (cosh ° ] cos(a;t). 
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We look for a sinusoidal solution of the form 

pl{t, W) = ^ [5(Vr)e'("*+<^(^))] , (24) 
with a real B{W) > 0. Inserting this form in Eq. fl23|) yields for (p and B 

W/ksT 

tan<p{W) = 1 (cos0<O) (25) 

2 cosh{So/ kBT) 



B{W) 



[cosh(5'o/A;BT)] 



-1 



[A cosh\So/ kBT) + (wro)2e2M^/'=s^]i/2- 

We now wish to quantify the hysteresis observed in the plane (/(t), U{t)). We choose a value 
Ui of the electric potential, such that < Ui < Uq. In the time interval —ir/u < t < tt/uj, 
there are two times, ti < and t2 = — ti such that U(ti) = U{t2) = Ui. Then the current 
intensity difference A/ = 1(^2) — I{ti) is a measure of the time variation of the conductance, 
since A/ = UiAG, with AG = ^(ts) - G{ti). From Eq. AG is given by 

AG = n,((72 - ^i)[pr(ti) - pr(t2)]. (26) 

To compute this last expression, we start by considering a single channel, that is, a fixed 
value of W. The difference Api{W) = Pi(ti, W^) — ^1(^2, W^) can be easily determined: 

Api{W) = ^^[cos{uh + 0) - cos(a;t2 + 4>)]- (27) 
ksT 

Taking into account the relation ^2 = — ^i, we get 

Ap,{W) = ^^M^sm<p{W)smut„ (28) 
kB-l- 

Evaluating sin0 from Eq. (1251) . we obtain 

^ 2B{W)au,ujr,e^l^^^{l-Ul/ulY/- 
' kBT[Acos\i\So/kBT) + {uToye^^l^BTY/2- V ) 

To obtain the current difference for the whole sample, we need to average over the energy 
barrier W: 

AI = Uin,{g2 - gi){Api)w. (30) 

With the notations 

h = 2 cosh{So/ kBT) (31) 
kBT cosh{So/kBT) \ J 
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we have 

(^P')- = " i, ''^•''^'^^ t^+(Jro)VW..T - (33) 

Introducing the change of variable x = cotq e^^''^'^ , we find 

'^''"-''LPi'""''^)^^ '''' 

with Tmin = roc^™/'^^"^. For cu <^ r^^, kBT\ii{x/ujTo) is typically large, and one can use 
the asymptotic expression f|T7|) of piW), yielding 

(35) 

If yU < 1, the integral converges to a finite value when its lower bound goes to zero, and 
we get that {Api)w scales as w^. The remaining integral can be computed exactly, and we 
eventually obtain for the average current variation 

^ 21+A'cosf [cosh(5o/A;BT)]2+M ^ ' 

xUiUo[l \\ n^{gi - g2) {ujtqY 

V '^0 J 

so that AJ also scales as u^, in the regime u <^ t^^ and /_i < 1 of experimental interest. 



D. Comparison between model and measurements 

Both calculations, in the case of a constant applied voltage and in the case of a varying 
voltage, lead to a power-law dependence of the conductance as a function of time in the 
first case and of frequency in the second. The model is therefore in qualitative agreement 
with our observations shown on Fig. 3 of a power law dependency. The experimental values 
obtained for the exponents m and m' should, according to the model, be equal, whereas 
they slightly differ -see Fig. 3-. Nevertheless, the large error bars on these experimental 
values make that they are still compatible with our model. 

Moreover, the model supposes that the exponent /i is equal to ksT /Wq, ie proportional 
to the temperature T. Indeed, we find a linear relation between the exponent m' and tem- 
perature (see Fig. 5). 

From our experimental observations we thus can roughly evaluate Wq, the "typical" value 
of the barrier height in the double-well model. Indeed, we have ii = ksT /Wq, so identifying 
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FIG. 5: (color online) The power law exponent m' as a function of T, obtained from fits of A/ as 
a function of oj. AI is evaluated at +400mV on a ±0.6^ I{U) cycle. A linear fit of the data is 
given. 

/i to m' gives Wq = 57meV. This value is on the order of magnitude of the Coulomb Blockade 
energy for one electron in a tunnel barrier—, which is yet another argument in favor of a 
microscopic origin for the conductance modifications in terms of trapping and untrapping 
of electrons on defects in the barrier. This would be the origin of the resistance switching 
observed in our tunnel junctions. 

We have to notice that the hypothesis made in the calculations, in the case of a periodic 
excitation, ie auQ <^ ksT, is not experimentally justified: auo = bAmeV if we assess a 
at 0.135eV/V and if we take uq = OAV. This value is larger than ksT in the studied 
temperature range, which means that the linearized expression of Eq. ( l23i) should be regarded 
as an approximation. 

We cannot extract more quantitative information from the comparison with experimental 
observations: our model does not predict the absolute value of the resistance relaxation with 
time, which would require for instance the knowledge of gi,g2 and tq. 

To be complete, we note that, on the one hand we observe telegraphic noise, thus 
associated to one defect, and on the other hand we model the junction in terms of a 
large distribution of defects, which could look contradictory. In fact telegraphic noise 
is observed at low temperature and low voltage, i.e., below 170mV. For these values, 
we can suppose within our model that state 1 is dominant. It means that the traps 
inside the MgO barrier are empty, all but one. This leads to the telegraphic noise. 
It is sometimes observed with several levels thus involving different defects, yielding 
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the addition of two telegraphic noise signals. On the contrary, at higher voltage or 
higher temperature, both populations -i.e. electrons on states 1 and 2- are present and 
the telegraphic noise is smeared out due to contributions of many defects. This explains 
that we observe a continuous relaxation of conductance with time, without telegraphic noise. 

V. CONCLUSION 

In conclusion, we showed here that a simple statistical model of electron trapping inside 
the MgO barrier could explain the resistance switching effects in MgO-based tunnel junc- 
tions. It also explains the long time relaxation of conductance according to a power-law 
behavior. In addition, the temperature dependence of the theoretical exponent is consistent 
with experimental observations. Our model supposes a change of the tunneling probability 
of electrons due to a local charging of the barrier: in that sense, it differs from usual hopping 
models through traps. 

We have to stress that this statistical model is in qualitative agreement with the phe- 
nomenological model that we proposed in Ref. [15]: in this model, inspired from memristor 
models^l, we introduced an electromigration term that makes the tunnel barrier height or 
thickness change as a function of the applied voltage, as well as an additional term which is 
voltage-independent. This second term makes the conductance relax towards a given value, 
independently of the applied voltage bias. Roughly speaking, this extra term plays a role 
analogous to that of thermal excitations in our present statistical model: thermal excita- 
tions indeed tend to equalize the populations of trapped and untrapped defects, and thus 
also tend to bring back the conductance to a given value, whatever the applied voltage. 

Our present approach is quite general for resistance switching effects in tunnel junctions: 
many resistance switching effects attributed for instance to ferroelectricity in the barrier 
-see for instance Ref. [27]- could perhaps be interpreted in terms of electron trapping on 
defects in the barrier. 
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